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With a sample of about 227 million BB pairs recorded with the BABAR detector at the PEP- 
II storage ring we perform a full angular analysis of the decay B° — > 0i^* (892) . We measure 
the branching fraction to be (9.2 ± 0.9 ± 0.5) x 10 -6 and determine the fractions of longitudinal 
and parity-odd transverse contributions as /l = 0.52 ± 0.05 ± 0.02 and f± — 0.22 ± 0.05 ± 0.02, 
respectively. The phases of the parity-even and parity-odd transverse amplitudes relative to the 
longitudinal amplitude are found to be <j> n = 2.34toil ± 0.05 rad and <j>x = 2.47 ± 0.25 ± 0.05 
rad, respectively. We measure five CP asymmetries which provide important limits on CP violation 
originating from new physics. We also observe the decay B° — > 4>K*°(1A30). 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The decay B — > tpK* (892) is expected to have contri- 
butions from b — > s loop transitions while the tree-level 
transition is suppressed in the Standard Model. Angular 
correlation measurements and asymmetries are particu- 
larly sensitive to amplitudes arising outside the Standard 
Model Q . The first evidence for this decay was provided 
by the CLEO and BABAR Q experiments. The large 
fraction of transverse polarization observed by BABAR |j| 
and confirmed by BELLE enables a full angular anal- 
ysis described by ten parameters for contributing ampli- 
tudes and their relative phases. 

The angular distribution of the B — > <f>K* decay prod- 
ucts can be expressed as a function of Hi = cos 9i and $, 
where 9i is the angle between the direction of the K from 
the K* — > Kit (6\) or (j> — > KK (0 2 ) and the direction op- 
posite the B in the vector resonance rest frame, and <& is 
the angle between the two resonance decay planes. The 
differential decay width has three complex amplitudes 
A\ corresponding to the vector meson helicity A = or 
±1 0, @. When the last two are expressed in terms of 
A|| = (A +1 + A_i)/V2 and A± = (A +1 - we 
have 

8tt d 3 r 1 

9T dHidH 2 d<S> ~ iAoPTp^FTpjTp X { 

\M 2 n\nl + \ (|A|,| 2 + \a^) (i - n\) (i - n\) 

+\(\A^ - \A^) [l -Hi) (l- Hi) co^ 
-lm{A ± Af) (1 - H\) (1 - Hi) sin 2$ 
+V2 ReiA^A*) Hi H 2 s] 1 - H\ sfl - H\ cos$ 
-V2Im(A±A* ) HxHi^l-H\^\-H\ sin* }. (1) 



In this analysis, we measure the branching fraction, 
obtained from the number of reconstructed signal events 
n s i g , the polarization parameters /l = |^4o| 2 /^I^a| 2 , 
f± = | 1 2 /S | ^4_a | 2 , and the relative phases 0|| = 
arg(A||/A ), 4>x = axg(Ax/A ). We allow for CP- 
violating differences between the B° (Q = +1) and B° 
(Q = —1) decay amplitudes (A\ and A\), where the fla- 
vor sign Q is determined in the self-tagging final state 



with a K* or K*: 

n% = n sig {l + QA C p)/2; (2) 
= h{l + QA° CP ); f± = f±Q- + Q*cp>, 
^ = </ }{] +QA^ ] ^ = ^ + | + Q(A^ + |). 

If one loop diagram dominates the decay amplitude, 
the three direct CP asymmetries Acp, A^ P , and Aq P , 
and the two weak-phase differences A0|| and A<f>±_ are 
expected to be negligible. From the above parameters 
one can derive vector triple-product asymmetries A\ and 
A'j* as discussed in Ref. 0: 

„||,o 1 ( MAxAj ) | Im(Z x I|, ) ^ 

T -*\ mx\ 2 + nz x \2 )■ [3) 

We use data collected with the BABAR detector Q at 
the PEP-II asymmetric-energy e + e~ collider || operated 
at the center-of-mass (CM) energy of the T(4S) reso- 
nance (y/s = 10.58 GeV). These data represent an in- 
tegrated luminosity of about 205 fb _1 , corresponding to 
226.6 ± 2.5 million BB pairs. 

Charged-particle momenta are measured in a tracking 
system consisting of a five-layer double-sided silicon ver- 
tex tracker and a 40-layer central drift chamber, both 
immersed in a 1.5-T solenoidal magnetic field. Charged- 
particle identification is provided by measurements of the 
energy loss (dE/dx) in the tracking devices and by a ring- 
imaging Cherenkov detector. 

We fully reconstruct B — > <fi K *° candidates from 
their decay products 6 — > K + K~ and K *° — > K^irT 
as discussed in Ref. |4j. Charged track candidates are 
required to originate from a single vertex near the inter- 
action point. We identify B meson candidates kinemat- 
ically using the beam-energy-substituted mass toes — 
[(s/2 + p t • p B ) 2 /Ef — pj] 1 / 2 and the energy difference 
AE = (EiE B - PrPs - s/2)/^/s, where (E h pi) is the 
initial state four-momentum obtained from the beam mo- 
menta, and (Eb,Pb) is the four-momentum of the re- 
constructed B candidate. The requirements on the K* 
and (j) invariant masses are 0.75 < toj^ < 1.05 and 
0.99 < m K x < 1.05 (GeV). We move the selection win- 
dow to 1.13 < niK-K < 1-73 (GeV) in the study of the 
higher-mass K* resonances. 
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To reject the dominant quark-antiquark continuum 
background, we require |cos$t| < 0.8, where 6t is the 
angle between the -B-candidate thrust axis and that of 
the rest of the tracks and neutral clusters in the event, 
calculated in the CM frame. We also construct a Fisher 
discriminant, T , further discriminating between signal 
and background, that combines the following variables: 
the polar angles of the B-momentum vector and the B- 
candidate thrust axis with respect to the beam axis in the 
CM frame, and the two Legendre moments Lq and £2 of 
the energy flow around the -B-candidate thrust axis [9|. 

Contamination from other B decays is small (about 
2% of the total background) according to Monte Carlo 
(MC) simulation and is taken into account in the fit 
described below. We remove signal candidates that have 
decay products with invariant mass within 12 MeV of the 
nominal mass values for Df 1 or — > (jm^ . 

We use an unbinncd, extended maximum-likelihood fit 
to extract simultaneously the signal yield and angular 
distributions from a sample of selected events. There are 
several event categories j: signal, continuum qq, combi- 
natoric BB background, B — > <fiKir with a non-resonant 



S-wave K^tt^ contribution, and B 
a broad S-wave K + K~ contribution 
each candidate i is defined as C 



* /o(980)A* with 
The likelihood for 

= E J , k ^Vf(x' t ;a-J), 

where each of the Vj(xi;a; (3) is the probability density 
function (PDF) for variables Xi = {wes, AS, T ', itikk, 
m KK> ^i) Q}- The flavor index k corresponds to 

the measured value of Q, that is Vj = Vj x 5kQ- The is 
the number of events with the flavor k in the category j. 

The PDF Vj(xi;a;P) for a given candidate i is the 
product of the PDFs for each of the variables and a joint 
PDF for the helicity angles and resonance masses as dis- 
cussed below. The signal angular distributions are pa- 
rameterized with the set a = {/l, f±, <j>i\, 4>±, A° CP , 
Aq P , A0||, A(f>±} which are left free to vary in the fit. 
The other PDF parameters j3 are extracted from MC 
simulation and data in jjjes and AE sidebands and are 
fixed in the fit. The MC resolutions are adjusted by com- 
paring data and MC in calibration channels with simi- 
lar kinematics and topology, such as B° — > D~ir + with 



D~ 



A" 4 



The PDF parameterization for each 



event candidate accounts for the loss of acceptance near 
Hi = 0.8 due to the Df and rejection requirements. 

We use a three-dimensional description for the helicity 
part of the signal PDF, using the ideal angular distri- 
bution from Eq. Q multiplied by an acceptance func- 
tion Q{H\, 7^2, 3>) parameterized with empirical polyno- 
mial functions. The detector acceptance effects arc found 
to be uniform in $, and we factor the Hi and Hi depen- 
dence as Q = Gi(Hi) x Qi^Hi). We use two Gaussian 
functions for the parameterization of the signal PDFs for 
AE, ttieS) and T . A relativistic P-wave Breit-Wigner 
distribution, convoluted with a Gaussian resolution func- 
tion, is used for the resonance masses. 



TABLE I: Summary of the 5° -> 0A*°(892) fit results. We 
show results for the ten primary signal fit parameters de- 
fined in Eq. @ and the derived parameters: reconstruction 
efficiency e which depends on decay polarization, branching 
fraction B, and triple-product asymmetries from Eq. Q. All 
results include systematic errors, which are quoted following 
the statistical errors. For the dominant correlations we give 
the coefficients in the last column. 



Fit parameter 



Fit result 



Correlation 



n a i g (events) 

h 
f± 

<t>\\ (rad) 
<t>± (rad) 
Acp 
A° CP 
Ab P 
A0|| (rad) 
A<j> ± (rad) 



201 ± 20 ± 6 
0.52 ±0.05 ±0.02 
0.22 ±0.05 ±0.02 

2 - M -o 2 20 ± °- 05 
2.47 ±0.25 ±0.05 
-0.01 ±0.09 ±0.02 
-0.06 ±0.10 ±0.01 
-0.10 ±0.24 ±0.05 
0.27^23 ±0.05 
0.36 ±0.25 ±0.05 



} -46% 
} ±70% 

} -45% 
} ±70% 



e(%) 
B 

A% 



9.7 ±0.5 
(9.2 ±0.9 ±0.5) x 10~ 6 
-0.02 ±0.04 ±0.01 
±0.11 ±0.05 ±0.01 



Parameterization of the non-resonant B-decay contri- 
butions is identical to that of the signal for ttt-es, AE, 
and T , but is different for the angular and invariant mass 
distributions. In particular, a broad invariant mass dis- 
tribution accounts for all potential S-wave contributions 
leaking into the mass selection window. For the com- 
binatorial background, we use polynomials, except for 
mES and T distributions which are parameterized by an 
empirical phase-space function and by the two Gaussian 
functions, respectively. Resonance production occurs in 
the background and this is taken into account in the 
PDF. The background Hi distribution is separated into 
contributions from combinatorial background and from 
real vector mesons. 

We allow for multiple candidates in a given event by 
assigning to each a weight of 1 /iVj , where iVj is the num- 
ber of candidates in the same event. The average number 
of candidates per event is 1.04. The extended likelihood 
for a sample of A can d candidates is 



£ = exp - n s Yl 



exp 



In C % 



(4) 



3 ! > 

The event yields nj , asymmetries Aj , and the signal po- 



m ES (GeV) AE (GeV) 




0.99 1.02 1.05 0.75 1.20 1.75 



m KK ( GeV ) m K7T (GeV) 

FIG. 1: Projections onto the variables rn.ES (a), AE (b), 
m KK ( c )i an d iriflSr (d) for the signal B° — > <f>K*° (892) and 
<^K"*°(1430) candidates combined. 




Q $ (degrees) Q 3> (degrees) 



FIG. 2: Projections onto the variables Hi (a), H2 (b), Q $ 
(c), and the differences between the Q<I> projections for events 
with Hi H 2 > and with Hi H 2 < (d) for the signal 5° -> 
<?!)if* (892) candidates. 



larization parameters a are obtained by maximizing C. 

The results of our maximum likelihood fit to the sam- 
ple of B° — > 0A"*°(892) candidates are summarized in 
Table [I] We also repeat the fit with the requirement 
1.13 < niKir < 1-73 (GeV) and without the angular in- 
formation. We observe 181 ± 17 events (statistical errors 
only) of the decays B° -> 0A*°(143O) with statistical 
significance greater than lOcr. In Fig.m^Jwe show pro- 
jections onto the variables, where data distributions are 
shown with a requirement on the signal-to-background 
probability ratio T-'sig/'Pbkg calculated with the plotted 
variable excluded. The solid (dashed) lines show the 
signal-plus-background (background) PDF projections. 

In the analysis of the decay B° -> </>A*°(892) 
for any given set of values (</>|| , cj>±, A(f>\\ , A<fi±) simple 
transformations of the angles, for example (— 4>\\,tt — 
4>±, —A<f)\\ , — A<j>±), give rise to other sets of values which 
satisfy Eq. Q in an identical manner. To resolve this am- 
biguity, the set of values lying closest to the theoretical 



FIG. 3: Projections onto the variables Hi (a) and H2 (b) for 
the signal B° -> A" (1430) candidates. The difference be- 
tween the solid and dotted lines in (a) shows the contribution 
of the tensor state to the angular distribution. 

expectation (tt, tt, 0, 0) 0,0, El * s cnosen - I n Fig. 0] we 
show likelihood function contour plots. 

We find the decay B° -> cf)K* a (U30) to be predomi- 
nantly longitudinally polarized based on the TL2 angular 
distribution in Fig. (b) . The width and the angu- 
lar distribution of the A*°(1430) resonance structure are 
not consistent with the pure A2°(1430) tensor state at 
more than 10a\ However, the angular distribution pro- 
vides evidence (with statistical significance of 3.2a) of the 
longitudinally polarized tensor (1430) contribution in 
addition to the scalar jFQf°(1430), see Fig. 0(a). 

Our B Q -> 0A*°(892) fit is performed with the B -> 
foK* and B — ► tpKn contributions unconstrained. We 
obtain the event yields 25 ± 10 and 11 ± 15, respectively. 
The systematic uncertainties due to interference are es- 
timated using generated samples with conservative as- 
sumptions about the S-wave intensity and the interfer- 
ence phase. Additional systematic uncertainty originat- 
ing from B background is taken as the difference between 
the fit results with the combinatoric BB background 
component fixed to zero and fixed to the expectation 
from MC. 

We vary the PDF parameters within their respective 
uncertainties, and derive the associated systematic er- 
rors. The biases from the finite resolution of the helicity 
angle measurement and the dilution due to the presence 
of fake combinations are estimated with MC simulation. 

The systematic errors in efficiencies are dominated by 
those in track finding and particle identification. Other 
systematic effects arise from event-selection criteria, 4> 
and A* branching fractions, MC statistics, and num- 
ber of B mesons. We calculate the efficiencies using the 
measured polarization and assign a systematic error cor- 
responding to the total polarization uncertainty. We find 
the uncertainty in the charge asymmetry due to the track 
reconstruction and identification to be less than 0.02. 

In summary, we have performed a full angular analysis 
and searched for CP violation in the angular distribution 

with the decays i?' -> < A > *°(892). Our results are 
summarized in Tabic [I] We observe, with more than 
5<7 significance, non-zero contributions from all of the 
three amplitudes \Ao\, \A±\, and \An\, see Fig.0|(a). We 
also find 3cr evidence for non-zero final-state-interaction 
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h 



FIG. 4: Likelihood function contours with la intervals for 
polarization (a) and phase (b) measurements in the B° — * 
<f>K*°(892) analysis. The fit results are shown with dots. Di- 
agonal dashed lines f± = (1 — /l)/2 and cj>± = (j>^ correspond 
to L4 + i 2> In (b) the (7r,7r) point is indicated by the 

crossed dashed lines. 



phases, see Fig. SI (b). These results supersede our earlier 
measurements in this channel 0, 0] . We also observe the 
decay B° -> 0Jf*°(143O). 

For B decays to light charmless particles we expect 
the hierarchy of decay amplitudes to be \A \ ^> \A +1 \ ^> 
\A-i\ under the assumption of pure loop diagram contri- 
bution, which is analogous to the discussion in Ref. 
Our measurements with the decay B° — > cj>K* Q (892) 
do not agree with the first inequality but agree with 
the previous measurements in Ref. 0, 0- This sug- 
gests other contributions to the decay amplitude, pre- 
viously neglected, either within or beyond the Standard 
Model [HQ. 

We would like to thank Yuval Grossman, Alexander 
Kagan, David London, Mahiko Suzuki, and Arkady Vain- 
shtein for useful discussions. We are grateful for the ex- 
cellent luminosity and machine conditions provided by 
our PEP-II colleagues, and for the substantial dedicated 
effort from the computing organizations that support 
BABAR. The collaborating institutions wish to thank 



SLAC for its support and kind hospitality. This work is 
supported by DOE and NSF (USA), NSERC (Canada), 
IHEP (China), CEA and CNRS-IN2P3 (France), BMBF 
and DFG (Germany), INFN (Italy), FOM (The Nether- 
lands), NFR (Norway), MIST (Russia), and PPARC 
(United Kingdom). Individuals have received support 
from CONACyT (Mexico), A. P. Sloan Foundation, Re- 
search Corporation, and Alexander von Humboldt Foun- 
dation. 



[1 

[2 
[3] 
[4] 
[5 

IT 

[8 

[9 

[10 

[11 

[12 
[13 



Now at Department of Physics, University of Warwick, 
Coventry, United Kingdom 

Also with Universita della Basilicata, Potenza, Italy 
Also with IFIC, Instituto de Ffsica Corpuscular, CSIC- 
Universidad de Valencia, Valencia, Spain 
Deceased 

G. Valencia, Phys. Rev. D 39, 3339 (1989); A. Datta and 
D. London, Int. J. Mod. Phys. A 19, 2505 (2004). 
CLEO Collaboration, RA. Briere et al, Phys. Rev. Lett. 

86, 3718 (2001). 

BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 

87, 151801 (2001); Phys. Rev. D 65, 051101 (2002). 
BABAR Collaboration, B. Aubert et al, Phys. Rev. Lett. 
91, 171802 (2003); arXiv:hep-ex/0303020 

BELLE Collaboration, K.-F. Chen et al, Phys. Rev. 
Lett. 91, 201801 (2003). 

G. Kramer and W.F. Palmer, Phys. Rev. D 45, 193 
(1992); H.-Y. Cheng and K.-C. Yang, Phys. Lett. B 511, 
40 (2001); C.-H. Chen, Y.-Y. Keum, and H-n. Li, Phys. 
Rev. D 66, 054013 (2002). 

BABAR Collaboration, B. Aubert et al, Nucl. Instrum. 
Methods A479, 1 (2002). 

PEP-II Conceptual Design Report, SLAC-R-418 (1993). 
BABAR Collaboration, B. Aubert et al, to ap pear in 
Phys. Rev. D (2004), arXiv:hep-ex/0403025 SLAC- 
PUB-10381. 

The BABAR detector Monte Carlo simulation is based on 
GEANT4: GEANT4 Collaboration, S. Agostinelli et al, 
Nucl. Instr. Meth. A 506, 250 (2003). 
A. Ali et al, Z. Phys. C 1, 269 (1979) 
Rev. D 66, 054018 (2002). 
S. Eidelman et al, Phys. Lett. B 592, 
Y. G rossman, Int. J. Mod. Ph ys. A 19, 907 (2004), A. Ka- 
gan, arXiv:hep-ph/0405134 P. Colangelo, F. De Fazio 
and T.N. Pham, arXiv:hep-ph/0406162 



M. Suzuki, Phys. 
1 (2004). 



